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Abstract. We investigate the supercluster MS0302+17 (z w 0.42) using weak lensing analysis and deep wide field 
BVR photometry with the CFH12K camera. Using (B — V) vs. (V — R) evolution tracks we identify early-type 
members of the supercluster, and foreground ellipticals. We derive a R band catalogue of background galaxies for 
weak lensing analysis. We compute the correlation functions of light and mass and their cross-correlation and test 
if light traces mass on supercluster, cluster and galaxy scales. 

We show that the data are consistent with this assertion. The (^-statistic applied in regions close to cluster 
centers and global correlation analyses over the whole field converge toward the simple relation M/L = 300 ± 
30^70 (M/L)q in the B band. This independently confirms the earlier results obtained by Kaiser et al. (1998). 
If we model dark matter halos around each early-type galaxy by a truncated isothermal sphere, we find that a 
linear relation M oc L still holds. In this case, the average halo truncation radius is s* < 200 h^g kpc close to 
clusters cores whereas it reaches a lower limit of ~ 300 /if kpc at the periphery. This change of s* as a function 
of radial distance may be interpreted as a result of tidal stripping of early type galaxies. Nevertheless the lack of 
information on the spatial distribution of late-type galaxies affects such conclusions concerning variations of s*. 
Though all the data at hand are clearly consistent with the assumption that mass is faithfully traced by light 
from early-type galaxies, we are not able to describe in detail the contribution of late type galaxies. We however 
found it to be small. Forthcoming wide surveys in UV, visible, and near infrared wavelengths will provide large 
enough samples to extend this analysis to late- type galaxies using photometric redshifts. 

Key words, observations: clusters of galaxies - general - cosmology:large-scale structure of Universe - gravitational 
lensing 
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K98), C11604+4304/C11604+4321 (Lubin ct al. 2000) or 
RXJ0848. 9+4452 (Rosati et al. 1999) show that supcrclus- 
ters of galaxies are genuine physical systems where grav- 
itational interactions between clusters of galaxies prevail. 
There is however no conclusive evidence that superclusters 
are gravitationally bound systems. 

Numerical simulations show that supercluster proper- 
ties may be best characterised by the shape and matter 
content of the filamentary structures connecting neighbor- 
ing groups and clusters. Unfortunately, the physical prop- 
erties of these filaments are still poorly known, though 
their existence at both low and very high redshift seems 
confirmed by a few optical and X-ray observations (Mollcr 
& Fynbo 2001; Durret et al. 2003). 

An alternative approach has been proposed by K98, 
G02 and Clowe et al. (1998) who used weak lensing 
analyses. While gravitational lensing has been used on a 
large sample of clusters of galaxies (See e.g. Mellier 1999; 
Bartclmann & Schneider 2001, and references therein), its 
use at supercluster scales is still in its infancy. It has been 
pioneered by K98 and Kaiser et al. (1999) who used V 
and I band data obtained at CFHT to probe the mat- 
ter distribution in the MS0302+17 supercluster. Similar 
analyses were done later by G02 for the A90 1/902 system 
using wide field images obtained with the WiFi instru- 
ment mounted at the MPG/ESO 2.2 telescope in La Silla 
Observatory. Wilson et al. (2001b) applied similar tech- 
niques to "empty" fields. In these papers, the projected 
mass density, as reconstructed from the distortion field 
of background galaxies, has been compared to the light 
distribution on several scales. All these studies conclude 
that there is a strong relation between the light of early 
type galaxies the and dark matter distributions as if all 
the mass was traced by early type galaxies. 

Over scales of a few megaparsecs, K98 found that 
the dark matter does not extend further than the light 
emission derived from the early-type galaxies sample. 
According to their study, there is almost no mass associ- 
ated with late type galaxies. In contrast, G02 argued that 
there is some, but its fraction varies from one cluster to 
another, leading to a lower and more scattered M/L than 
K98 found. This discrepancy may be explained if either 
the two superclusters are at different dynamical stages, or 
their galaxy populations differ (fraction of early/late type 
galaxies) . 

The reliability of their results may however strongly 
depend on systematic residuals. It is worth noticing that 
the weak lensing signal produced by filaments is indeed 
expected to be difficult to detect because projection ef- 
fects may seriously dilute the lensing signal (Jain, Seljak, 
& White 2000). Therefore, systematics produced by tech- 
nical problems related to the way both groups analyzed 
their data may also be a strong limitation. It is therefore 
important to confirm early K98 conclusions from an inde- 
pendent analysis, and possibly go further in order to take 
into account the properties of dark halos. In particular, it 
is interesting to compare the halo properties (size, velocity 
dispersion) of cluster galaxies with those of field galaxies. 



MS0302+17 seems a generic and almost ideal z ~ 0.42 
supercluster configuration for such an astrophysical study 
because it is composed of three very close rich clusters 
(mean projected distance ~ 15'). 

In this paper, we describe the investigation of this sys- 
tem. Using new data sets obtained in B, V and R at 
the CFHT with the CFH12K CCD camera, we explore 
the mass and light distributions in the supercluster area. 
Since the CFH12K field of view is 1.7 times larger than 
the UH8k camera, the three clusters and most of their 
periphery are totally encompassed in the CFH12K field 
and we can even explore whether other clusters lie in the 
field at the same redshift. A quick visual inspection of the 
southern cluster (CIS) reveals that it is very dense. The 
giant arc discovered by Mathez et al. (1992) is visible. 
Similar arc(let) features are also visible in the northern 
(C1N) and eastern (C1E) systems, making the MS0302+17 
supercluster a unique spectacular lensing configuration, 
where strong and weak lensing inversions can be done. 

This paper is organized as follows. In Sect. 2 we de- 
scribe the observations that were carried out and re- 
duced. We include details on how both astrometric and 
photometric solutions were computed. We present a de- 
tailed quality assessment of the catalogues, where com- 
parisons to existing deep catalogues are made. In Sect. 3 
we explain how object catalogues were produced for su- 
percluster members identification. Sect. 4 presents the 
weak lensing signal produced by the dark matter com- 
ponent. Sect. 5 shows how these two components cross- 
correlate. Results are discussed in Sect. 6 and our conclu- 
sions and summary are presented in Sect. 7. Throughout 
this paper we adopt the cosmology fi m = 0.3, Ha = 
0.7, Hq = 70 /170 km.s _1 Mpc _1 , leading to the scaling 
relation 1' = 333 h^ kpc at z — 0.42. 

2. The Data 

2.1. Observations and data reduction 

The observations of the MS0302+17 supercluster area 
were obtained on October 12, 1999. They were carried 
out with the CFH12K camera mounted at the prime fo- 
cus of the Canada-France-Hawaii Telescope. The CFH12K 
mosaic device is composed of 6 x 2 thinned backside illu- 
minated MIT Lincoln Laboratories CCID20 2048 x 4096 
CCDs with a 15 /im pixel size. The wide field corrector in- 
stalled at the prime focus provides an average pixel scale 
of 0'.'205 and the whole field of view of the CFH12K cam- 
era is 43.'2 x 28.'9. Useful details on the camera can be 
found in Cuillandrc et al. (2000) and in McCracken et al. 
(2003) (hereafter McC03). 

The pointings were centered at the reference posi- 
tion RA(2000) = 03 h 05 m 26.00 s and DEC(2000) = 
+17° 17' 54", so that the CFH12K field of view encom- 
passes the three major clusters of galaxies. Sequences of 
dithered exposures were obtained in B, V and R filters us- 
ing small shifts of about 30 arcsec to fill the gaps between 
the CCDs and to accurately flat field each individual im- 
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Table 1. Summary of the observations with total expo- 
sure time (+ the number of dithered pointings), seeing and 
limiting magnitude. Following McCracken et al. (2003), 
the AB limiting magnitude corresponds to a 50% com- 
pleteness limit. 



Filter 


exp. time 


Seeing 


Limiting 






(seconds) 


(arcsec) 


Mag 




Bab 


6 x 600 = 3600 


0.9 


25.75 




Vab 


6 x 600 = 3600 


0.9 


25.50 


d 


Rab 


10 x 720 = 7200 


0.7 


26.50 


oc 

0) 



age frame. Table 1 summarizes and assesses the useful 
observations used in this paper. 

The B, V and R data were processed and calibrated 
at the TERAPIX data center located at IAP. The pre- 
calibration process was done using the flips package 1 . 
Photometric and astrometric calibrations, as well as im- 
age stacking and catalogue production were done using 
the current software package available at TERAPIX 2 . The 
overall pre-reduction (bias and dark subtraction, flat-field 
calibration), photometric and astrometric calibration as 
well as image resampling and co-addition follow exactly 
the same algorithms and steps as in McC03. We refer to 
this paper for further details. 

2.2. Photometric calibration 

The photometric calibration was done using Landolt 
star fields SA95 and SA113 (Landolt 1992). The IRAS 
maps (Schlegel, Finkbeiner & Davis 1998) show that the 
Galactic extinction is important in this field. The av- 
erage E(B — V) = 0.125, leading to extinction correc- 
tions in B, V and R of 0.508, 0.384 and 0.285 respec- 
tively We adopted the AB magnitude corrections pro- 
vided by McC03: B AB =B- 0.097, Vab = V- 0.007 and 
Rab = R + 0.218. 

The object photometry was derived using the Magauto 
parameter of SExtrator (Bertin & Arnouts 1996) for the 
magnitude and Magaper for the color index, with the V 
catalogue as reference position, inside a 2 arcsec aperture. 
The galaxy number counts derived from this photometry 
peak at B AB = 24.9, Vab = 24.3 and Rab = 25.6. As in 
McC03, we derived the limiting magnitude by adding sim- 
ulated stellar sources in the field. The limiting depths at 
which 50% of these sources are recovered are Bab = 24.9, 
Vab = 24.3 and Rab = 25.6, in good agreement with the 
expectations from the McC03 F02 deep exposures, once 
rescaled to a similar exposure time (see Fig. 1). 

The reliability of the photometric calibration was 
checked by comparing the B, V and R galaxy counts to 
those published for the VIRMOS F02 field by McC03. The 
supercluster galaxy population contaminates the galaxy 
number counts. However, as the B — (B — V) and R — 
(B — R) color-magnitude diagrams show, the bright-end 




1 http : //www. cf ht .hawaii . edu/~ jcc/Flips/f lips .html 

2 http://terapix.iap.fr/ 



AB mag. 

Fig. 1. BVR galaxy counts in the MS0302+17 superclus- 
ter field. The red, green and blue dots show the galaxy 
number density in R, V and JB, respectively. The solid red, 
green and blue lines show the same counts obtained with 
the same filters on the deep VIRMOS F02 field by McC03. 
For the bright and faint magnitudes the agreements are 
excellent. In between, the MS0302+17 counts show a sys- 
tematic excess that reveals the supercluster populations 
at z = 0.42. 

and faint-end galaxy populations are dominated by field 
galaxies, so we expect the data of the MS0302+17 su- 
percluster and the VIRMOS F02 data to be comparable 
for these populations. We checked that the counts agree 
within 0.05 magnitudes with McC03 for the three filters. 

We also checked the (B — V) versus (V — R) color-color 
magnitude of stars, as selected from the SExtractor stel- 
lar index. To avoid mixing of galaxies in the sample, we 
only used bright objects that have a reliable stellar index. 
The colors are plotted in Fig. 2 (yellow dots) and com- 
pared to selected stars on the main sequence and giant-star 
(B — V)/{V — R) tracks from Johnson (1966). Two popu- 
lations are clearly visible: the blue stars [B— V < 1.0) are 
halo stars, the red ones ((B — V) ~ 1.3) are disk M-dwarf 
stars. The blue population perfectly matches the Main 
Sequence, but for the red sample we found a systematic 
offset of 0.1 in the (B — V) term. A similar discrepancy 
has already been mentioned by Prandoni et al. (1999). It 
is likely due to the color correction, which is no longer 
linear for those stars. The dispersion is easily explained 
by the intrinsic dispersion of stellar colors and likely from 
statistical magnitude measurement errors. We therefore 
conclude that our (B — V) and (V — R) colors have an 
internal error of ±0.05 magnitude. 

The CCD to CCD calibration errors are in principle 
minimized since all CCDs are rescaled with respect to 
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Fig. 2. Color-color diagram of bright stellar objects in the 
MS0302+17 supercluster field. The horizontal and verti- 
cal axes are the Johnson (B — V) and (V — R) color in- 
dices, respectively. The red and blue dots show the stellar 
tracks for Main Sequence and Giant stars, as predicted by 
Johnson (1966). The yellow dots are the color positions of 
the bright star sample. The figure clearly shows two pop- 
ulations : the blue halo stars ((B — V) < 1.0) and the red 
disk M-dwarf stars. 



the reference CCD#4 that contains several Landolt stars. 
However, residuals from illumination correction may still 
bias the calibration. We checked the CCD to CCD sta- 
bility by comparing the B, V and R galaxy counts. The 
contamination by cluster galaxies in the field makes this 
approach difficult. Their effects were minimized by remov- 
ing the cluster regions from the count estimates. However, 
since the diffuse supercluster filaments also contaminate 
the signal, large CCD to CCD fluctuations of the counts 
still remain. We therefore focused on the faint end part of 
the magnitude distribution, where the supercluster popu- 
lations should be negligible. Using these constraints, the 
average CCD to CCD galaxy count fluctuations in B, V 
and R are 2.5% in each filter, with peaks of 7.5%. When 
clusters fields are included and the magnitude range is 
broadened, the peaks reach 16%. This clearly reveals the 
presence of clusters populations. Possible residuals from 
calibration problems are therefore negligible compared to 
fluctuations expected from Poisson or cosmic variance. 

Possible color variations of stars from one CCD to an- 
other were checked using a color analysis of 100 stars per 
chip. We found that they all show similar average (B — V) 
and (V-R) colors inside the 12 CCDs, within ±0.1 (iela) 
and without any systematic color gradient. This confirms, 
independently of the galaxy counts, that the photometry 



and our [B — V) and (V ~ R) colors are stable enough 
across the field to meet our scientific goals. 
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Fig. 3. Astrometric calibration of the MS0302+17 field. 
The top panel shows the residual Aa and Ai5 between 
the USNO reference star position and the astrometric so- 
lution. The two bottom panels show the RA and DEC 
residual star position difference as a function of the RA 
and DEC in the MS0302+17 supercluster field. No trend 
is visible, showing that astrometry is stable across the 
whole field of view. 



2.3. Astrometric calibration 

The astrometric calibration was done using the 
Astrometrix package developed jointly by the TERAPIX 
center and Osservatorio di Capodimonte in Naples for 



Gavazzi et al.: Mass and Light in the MS0302+17 Supercluster 



5 



wide field images 3 . The algorithms are extensively de- 
scribed in McC03 and Radovich ct al. (2003), so we refer 
to these two papers for further details. 

For the MS0302+17 data, the calibration was done us- 
ing the USNO-A2.0 reference star catalogue (Monet 1998). 
The astrometric center and tangent point of the CFH12K 
is a 200 o = 03 h 05 m 25.8 s and <5 2 ooo = +17° 17' 54". Since 
the V image has much fewer saturated stars than the R 
image, we preferred to use it as the reference astrometric 
data, although the R is deeper and has a better seeing. 
We then cross-correlated the V input catalogue generated 
by SExtractor to the USNO-A2.0. After rejection of sat- 
urated objects, we found 731 stars common to both cat- 
alogues. The star sample is homogeneously spread over 
the CCDs, so that we used 50 to 70 stars per CCD. Fig. 3 
shows the residuals between the reference USNO-A2.0 star 
positions and positions derived from the astrometric solu- 
tion. The rms coordinate error is 0'.'45 (68%) and is similar 
in both RA and DEC directions (Aa = 0'.'36, A<5 = 0'.'30, 
respectively). No systematic shift or position gradient is 
visible. Each V frame was resampled according to the as- 
trometric calibration and then stacked to produce the final 

V image. 

The B and R images are calibrated with respect to the 

V image. The matching uses detection catalogues gener- 
ated by SExtractor so the cross-correlation can be done 
using several thousands of stars and galaxies. Since the 
B, V and R data were obtained during the same night, 
systematic offsets only correspond to small shifts imposed 
by the observer, and rotation between each image is neg- 
ligible. Across such a small field, the atmospheric differ- 
ential refraction for an airmass ranging between 1.0 to 1.5 
produces shifts smaller than 0'/4 between the B, V and 
R image and a chromatic residual that is less than 0'.'04 
between the B and R image. The rms coordinate errors 
between the B and R catalogue and the V reference are 
0'.'05 (68% CL). 

2.4. Final catalogues 

Using the calibrated B, V, and R images, we then pro- 
duced the final BVR objects catalogue. It contains 28600 
galaxies and 1100 stars, as defined according to the 
SExtractor stellarity index. The common area is com- 
posed of 12500 x 8500 pixels and is 0.343 deg 2 wide. 
It corresponds to an angular size of 42 '.7 x 28'.9 (i.e. 
14.2 Mpc x 9.6 Mpc at z =0.42) . 

For weak lensing analysis we produced another cata- 
logue that only uses the R band image without regard of 
the B and V data. This image is the deepest one and has 
the best seeing. The weak lensing catalogue contains more 
objects than the joined B, V and R one. Its properties are 
detailed in Sect. 4. 

CCD defects, gaps or overlap areas between CCDs, 
bright stellar halos, saturated stars and asteroid track 
residuals generate spurious features. They are removed 



from the catalogue using the manual masking process de- 
scribed in van Waerbeke et al. (2000, 2001) as shown in 
Fig. 4. Field boundaries are also masked, and we finally 
end up with an effective area of 0.228 deg 2 . The catalogues 
discussed in the following will only concern this common 
unmasked part of the field. 

3. Supercluster galaxies and light distribution 

As shown above, the color information is stable across the 
field within 0.05 mag accuracy which is sufficient to make a 
reliable selection of cluster and non-cluster galaxies using 
colors. The supercluster member selection and the redshift 
distribution of foreground and background lensed sources 
were done using color-color diagrams together with the 
measurement of photometric redshifts. 

3.1. Photometric Redshifts vs. Color-Color relation 

We first attempted to use B, V and R photometry to 
derive photometric redshifts using the hyperz 4 package 
(Bolzonella et al. 2000). Details of the method applied 
to clusters or deep multi-color wide field surveys can be 
found in Athreya et al. (2002); van Waerbeke et al. (2002); 
Gavazzi et al. (2003). However, compared to these previ- 
ous analyses, we only have three bands, which severely re- 
duces the reliability of photometric redshift information. 
When compared to an "empty" region located westward 
on the field, the photometric redshifts show an excess of 
galaxies in the redshift range [0.25 — 0.65]. The photomet- 
ric redshift uncertainty (~ ±0.2) hampers any detailed 
redshift investigation of the supercluster galaxies. This 
noisy redshift information can partially be used for the 
distinction between foreground and supercluster objects 
and background lensed galaxies (See Sect. 4). 

Fortunately, at redshift z w 0.4 the typical 
4000A break spectral feature lies between the B and 
R filters and can easily reveal early type galaxies. The 
cluster selection was therefore primarily focused on the 
red cluster sequence in the color-magnitude diagrams. By 
using a B, V and R color-color diagram of well-defined 
magnitude limited sample of galaxies for the selection the 
supercluster members, one can easily isolate co-eval early- 
type cluster galaxies. We first flagged objects within a 1.5 
arcmin radius from the center of each cluster. Fig. 5 shows 
the (B — V) versus (V — R) color-color diagram. Early- 
type galaxies concentrate around 1.1 < (B — V) < 1.42, 
the supercluster elliptical galaxies at z = 0.42 having also 
0.47 < {V - R) - 0.5(5 - V) < 0.77. Objects with a 
(V — R) color bluer than that band are elliptical galax- 
ies with a lower redshift. Most are in the redshift range 
0.1 < z < 0.4. Passive evolution tracks for elliptical galax- 
ies kindly provided by D. Leborgne (see Fioc & Rocca- 
Volmerange 1997) are in excellent agreement with our ob- 
servations. 



3 http : //www . na . astro . it /"radovich/ 



4 http : / / webast . ast . obs-mip . f r/hyperz/ 
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Fig. 4. i? band image of the MS0302+17 supercluster field used for the light and mass analysis. Coordinates are in 
J2000. The center point is RA = 03 h 05 m 25.79 s and DEC = +17° 17' 54.02". The field size is 42.9' x 29.0' (i.e. 
14.2 Mpc x 9.6 /ifo 1 Mpc @ z = 0.42). The 3 arcmin radius circles show the cluster positions of the X-ray emission 
peaks: C1N, North 03 h 05 m 17.81 s , +17° 28' 37.6", z = 0.425, C1E, East 03 h 06 m 18.89 s , +17° 18' 33.9", z = 0.418 
and CIS, South 03 h 05 m 31.49 s , +17° 10' 16.3", z = 0.426. Green polygons define the masked areas that have been 
removed from the image analysis. 



The following set of equations summarizes our selec- 
tion criteria for early-type galaxies. Equations (la)-(lc) 
are criteria for supercluster members, whereas (la)-(lb) 
and (Id) stand for the selection of foreground early- type 
galaxies. 

19 < R < 23, (la) 
1.1 < (B-V) < 1.42 (lb) 
0.47 < (V - R) - Q.5(-B - V) < 0.77 (lc) 
0.2 < (V-R)- 0.5(5 - V) < 0.47 (Id) 

The foreground sample has a mean redshift z ~ 0.3 
and contains 770 galaxies, and the supercluster sample 
contains 750 galaxies. Their luminosity distribution is well 
centered around L*. We did not select ellipticals with 
z < 0.2 since they are very few and provide a negligi- 
ble lensing signal. Elliptical galaxies at z > 0.6 also have 
a poor lensing efficiency compared to those at the super- 
cluster redshift. 

Our color-color selection method fails to localize bluer 
late-type galaxies. Thus, we have to keep in mind that 



the light due to cluster spirals is not taken into account. 
K98 and Fabricant et al. (1994) argued that « 30% of the 
total B band rest-frame luminosity is due to late-types, so 
their contribution, though sizeable, is not expected to be 
dominant. 



We estimated the contamination by field galaxies in- 
side the color-color region of supercluster members. We 
selected similar galaxies satisfying (la)-(lc) far from the 
regions around the three clusters and plotted them in the 
color-color diagram. The fraction of field galaxies inside 
the cluster color-color region turns out to be negligible. 
This efficient selection process expresses the fact that for 
this redshift (z ps 0.4), the (B - V) and (V - R) colors are 
reliable filters. In contrast, the lower (V—R) limit used for 
the foreground subsample selection is more questionable, 
so we may miss some of the nearest ellipticals. However, 
since their lensing contribution is small, it has no impact 
on the interpretation of the weak lensing signal. 
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Fig. 5. Upper left panel: Color-color diagram. The red + 
signs code for the supercluster objects that lie within 1.5 
arcmin of individual clusters. The early-types are concen- 
trated around (B-V)ab ~ 1-2. The additional knowledge 
of (V — R) color provides an efficient selection of super- 
cluster objects at z ~ 0.4. Upper right panel: Color-color 
diagram that explicitly makes the distinction between su- 
percluster and foreground ellipticals along the (V — R) 
axis. The variation of (V — R) with redshift along the 
(B — V)ab ~ 1.2 line is well reproduced by predicted pas- 
sive evolution tracks for ellipticals whatever their forma- 
tion redshift. We define a foreground subsample centered 
around z w 0.3 i.e. inside the lower green lozenge whereas 
the upper red one encompasses the supercluster ellipticals. 
The lower panel shows the corresponding spatial distribu- 
tion. Supercluster members (red -I- symbols) are clearly 
clustered around the three known clusters. Foreground ob- 
jects (green x symbols) are clearly less clustered although 
some clumps are also visible. 



3.2. Spatial distribution of light & associated 
convergence 

In this section, we investigate the spatial distribution of 
supercluster galaxies. We also estimate the foreground 
contribution to light. The apparent Bab magnitudes are 
converted into rest-frame luminosities using K-corrections 
also provided by D. Leborgne. For z = 0.42 superclus- 
ter members, we used Kb = 1-88, and for z ~ 0.3 fore- 



ground early-types Kb — 1.14. As detailed in Sect. 4.2, 
the redshift distribution of lensed sources implies the 
following critical surface densities: Yj cr i t {zd — 0.42) = 
2.72 x 10 14 M Q /arcmin 2 and £ crit (z d = 0.3) = 1.95 x 
10 14 h^M® /arcmin 2 . 

Instead of computing a filtered luminosity map with 
Gaussian smoothing, we follow the method of K98, G02 
and Wilson et al. (2001b). It consists of weighting the 
luminosity map by its lensing efficiency (different for fore- 
ground and supercluster components). Assuming a con- 
stant mass-to- light ratio for supercluster galaxies, the lu- 
minosity field at a given position 6 is converted into a mass 
density field, Y.c(9). Then it is translated into a conver- 
gence Kc(0) = S£(0)/E cr it field. We finally convert it into 
a shear field j{9) and derive a supercluster shear pattern 
that samples the field according to the source catalogue 
positions (see Sect. 4.2). From this shear field a new con- 
vergence field can be drawn. It has the same field size and 
shape, the same masking and the same sampling proper- 
ties as the K-map we construct in Sect. 4.3. 

The inferred shear field reads : 



7 (0) 



£70(0 



9i) 



W{9), 



(2) 



where IF is a 40 arcsec Gaussian smoothing filter, and 70 
is the shear profile of an individual galaxy. Since we are 
primarily interested in the collective behavior of galaxies 
and we are looking at scales > 1 arcmin, we do not make 
a further hypothesis about the radial density profiles of 
galaxy halos. This means that galaxies are equivalent to 
point masses, so 70 simply reads: 



7o (0) 



M 



L ttE 



crit c 



(3) 



where M/L is the mass-to-light ratio which is assumed to 
be the same for all galaxies. The 7 — ► k inversion is de- 
tailed in Sect. 4.3. It turns out that luminosity-weighted 
«£ and number-density- weighted 5 kn convergence maps 
are almost proportional. A small deviation from equal- 
ity appears for the highest contrast values. In this case, 
Akjv/ktv Ak£/k£ as we expect if the brightest galaxies 
lie in the densest regions. The resulting n from light maps 
are shown in Fig. 6. The upper panel shows the kc con- 
vergence map for the supercluster objects only, whereas 
the lower panel shows the modifications produced by ad- 
dition of foreground structures. The three known clusters 
are clearly detected and seem to encircle a large under- 
dense region. A diffuse extension, less dense than the clus- 
ters, appears westward from the northern cluster (C1N). 
This extension encompasses two clumps at Aa ~ — 8' and 
Aa ~ —15'. Another extension toward the North- West 
of CIS is partly due to foreground structures. Using the 
spectroscopic redshift of two member galaxies, K98 argued 
that this clump probably lies at z = 0.3 . 



in this case, 70 (9) 



Mo 



v gz with Mo the mass of a galaxy 
halo, which is assumed to be constant from one galaxy to an- 
other. 
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that bias cllipticity measurements. The reliability of shape 
measurements is expected to be as good as the current 
cosmic shear survey data (van Waerbeke et al. 2001). 

4.1. PS F correction 

Blurring and distortion of stars and galaxies produced 
by instrument defects, optical aberrations, telescope guid- 
ing, atmospheric seeing and differential refraction are cor- 
rected using the PSF of stars over the whole field. Several 
correction techniques and control of systematic errors have 
been proposed over the past 10 years (see e.g. Mellier 
1999; Bartelmann & Schneider 2001; van Waerbeke & 
Mellier 2003; Refregier 2003). In the following we use the 
most popular KSB95 method initially proposed by Kaiser, 
Squires, & Broadhurst (1995). Several teams have already 
demonstrated that the KSB95 method can correct system- 
atics residuals down to the lower limit shear amplitude 
expected on supercluster scales (van Waerbeke & Mellier 
2003; Refregier 2003). 

Following KSB95 method, the observed ellipticity com- 
ponents e°ifj 2 are composed of its intrinsic ellipticity 
components e^ rc , and linear distortion terms that express 
the instrument and atmospheric contaminations and the 
contribution of gravitational shear to the galaxy elliptic- 
ity. Each ellipticity component is transformed as: 



Fig. 6. Effective convergence maps derived from the lu- 
minosities of early type galaxies. The convergence is ex- 
pressed as k(0) — £(0)/£ crit but the supercluster con- 
tribution may be contaminated by foreground galaxies. 
Panel a) (top) shows the luminosity weighted conver- 
gence map of supercluster objects, as inferred from the 
positions of galaxies in the color-color diagram. Panel b) 
(bottom) shows the same map with the contribution of 
the foreground galaxies. The origin coordinate is RA = 
03 h 05 m 25.79 s and DEC = +17° 17' 54.02". We applied a 
40" wide Gaussian smoothing scale. The two maps are al- 
most similar, confirming that foreground structures along 
the line of sight do not dominate uncertainties in the error 
budget. Green levels start at k — 0.0 and increase linearly 
with step 0.02. We assumed a fiducial mass-to-light ratio 
M/L = 300/i70 (M/L) . 



4. Weakly Lensed objects sample - Shear analysis 

The coherent stretching produced by the weak lensing ef- 
fect due to the MS0203+17 supercluster is measured using 
the deep catalogue extracted from the R band image. Its 
depth and high image quality allow us to lower the detec- 
tion threshold and to increase the galaxy number density 
(ss 25 arcmin~ 2 ), compared to the B and V images. This 
reduces the Poisson noise of the weak lensing statistics and 
increases the spatial sampling of the supercluster mass re- 
constructions. Close galaxy pairs with angular separations 
less than 5 arcsec are discarded to avoid blended systems 



With i*. 



psh 



psm 

Q7 



psm * 
psh 



(4a) 
(4b) 



7/3 



where g is the reduced gravitational shear, P sm is the 
smear polarizability, P sh the shear polarizability and 
P s the isotropic circularization contribution to the final 
smearing. In the following, all these tensors are simplified 
to half their trace and have been calculated with Imcat 6 
tools, fjssrj and q* are quantities that are measured 
from field stars. Their shape is fitted by a second or- 
der polynomial, applied individually to each CCD of the 
CFH12K camera. Stars are selected in the magnitude-r^ 
plane, as usual, q* is the anisotropic part of the PSF, 
which is subtracted from galaxy ellipticities. The residual 
is shown in Fig. 7 . It does not show any peculiar spatial 
pattern and is consistent with a one percent rms noise. 

The smearing part of the PSF contained in the P s 
term depends on the magnitude of the object and on its 
size as compared to the seeing disk. To optimally extract 
P g , we derived it from an averaged value over its 70 near- 
est neighbors in the magnitude— plane. Its variance is 
then used as a weighting scheme for the shear analysis. 
The weight assigned to each galaxy is finally the inverse 
variance of. of the observed ellipticities : 



1 



ps2 



pg2 c 



(5) 



http : //www. if a.hawaii . edu/ "kaiser/ imcat/ 
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Fig. 7. Upper panel : Spatial distribution of stellar el- 
lipticities before (black) and after (red) PSF correction. 
Lower panel : The same stellar ellipticities in the (ei,e2) 
plane, uncorrected (black crosses) and corrected (red dots) 
from PSF anisotropy. The residual rms dispersion around 
the center is less than one percent. 



where ctq = 0.3 is the intrinsic dispersion in galaxy ellip- 
ticities and <7j « 0.13 is the observed dispersion of elliptic- 
ities over the 70 closest neighbors. A weighted magnitude 
histogram of the source sample shows that we can select 



galaxies down to the limiting magnitude Rab = 25.6. In 
this subsample, galaxies with a smaller than that of 
stars are discarded. 



4.2. Redshift distribution & shear calibration 

Because the source sample is deeper than the photometric 
catalogue, the redshift distribution of background objects 
cannot be derived from the photometric redshifts calcu- 
lated with the B, V and R data. Nevertheless, for the 
distinction between foreground and background-lensed ob- 
jects, we can use a rough limiting redshift estimate. 

We chose to reject objects with a 2 p hot < 0.5. We 
also discarded the early-type galaxies selected in the color- 
color diagram of Sect. 3. Finally, we selected background 
galaxies within the magnitude range : 22.4 < Rab < 25.6. 
The source catalogue contains 22125 galaxies. This corre- 
sponds to a number density rib g ~ 27arcmin -2 . 

The properties of the resulting sample are roughly 
comparable with those of the sample of van Waerbeke et 
al. (2001, 2002), though their magnitude cut I AB < 24.5 
instead of Rab ^ 25.6 as in this work. They inferred the 
redshift distribution: 



n(z) 



1 



z,rfa) 



— e 



-z/z s 



(6) 



where a = 5/2, z s = 0.44 leading to z = az s — 1.1 and 
a z = s/aZg = 0.7. The median redshift is well approx- 
imated by Med(z) ~ (a — 0.33)z s ~ 0.95. At the same 
time, G02 proposed a median redshift 2 = 1 for their 
magnitude cut R < 26. Using the same analytic form as 
(6), we found that a ~ 1.9 and z s ~ 0.55 provide a good 
description of our redshift distribution implying a median 
redshift Med(z) = 1 and a broader distribution <r z ~ 0.8 . 
For the supercluster redshift Zd = 0.42 we calculated the 
mean of the ratio [3 — (Dd s / D s ) and the corresponding 

2 a- 1 

critical surface density S cr jt = 4^75/57 1 where Dd, D s 
and Dd s are angular distances between the observer and 
deflector, observer and sources and deflector and sources, 
respectively. We found : 



f3 = 0.49, 
S crit = 2.72 x 10 15 h 70 M Mpc~ 2 

= 3.02 x 10 14 M arcmkr 2 . 



(7) 



For foreground galaxies at z ~ 0.3, we found (3 = 0.63, 
Sent = 1.95 x 10 14 M arcmin~ 2 . The redshift distri- 
bution of sources is indeed equivalent to a single source 
plane configuration with redshift z s hcet ~ 0.95. The depth 
and the source plane redshift we use are in good agree- 
ment with previous ground based analyses like those of 
Clowe & Schneider (2001, 2002). The uncertainty in the 
gravitational convergence produced by the redshift distri- 
bution of the sources is about 5%, which is much smaller 
than the error bars we expect from statistical noise due to 
intrinsic galaxy ellipticities. 
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4.3. Mass Map 

Our mass reconstruction is based on the Kaiser & 
Squires (1993) (KS93) algorithm. The convergence k{6) = 
£(0)/£ C rit is related to the observed shear field 7(0) 
through: 

k{0)= I L>(0-tf)*7(tf)d 2 tf, (8) 

where 7 and D{6) = ^ m^E^p are complex quantities. 
On the physical scales we are exploring the lensing signal 
is weak enough so that (e) = ~ 7. The ellipticity 
catalogue is smoothed with a 9 S = 40" Gaussian filter : 

TW^E^expf-^^V (9) 

where wt are the weights defined in Eq. (5) and N ~ 
27rnbg^g ~ 170 can be viewed as the mean number of 
sources inside the filter. The resulting convergence map 
presents correlated noise properties : 

^^^ )>^ eXP (- i lf)- (10) 

<7, sa 0.42 is the dispersion in ellipticities of our galaxy 
sample. , CT ' ~ 0.016 characterizes the noise level. 

-y/87rn bg 02 

The K-map reconstruction result is shown in the middle 
panel of Fig. 8. The bottom panel shows the reconstruc- 
tion applied to the same galaxy sample, but with the ori- 
entation rotated by 45°. It represents the imaginary part 
of Eq. (8), which should be a pure noise realization if the 
coherent distortion field is only produced by gravitational 
lensing. 

The three main clusters C1N, CIS & C1E are de- 
tected with a high significance. A few substructures with 
a lower detection significance are also visible. We detail in 
Sect. 4.4 the various quantities measured for the clusters. 
A comparison with the X-ray emissivity map displayed 
on Fig. 1 of K98 shows an excellent agreement. See also 
Fig. 1 of Fabricant et al. (1994). The clusters are also de- 
tected in the K-from-light map shown in the top panel of 
Fig. 8. The clumpy extension detected westward from C1N 
can be seen in the light map. Another clump located at 
03 h 04 m 30 s , +17° 15' 00" may be either a foreground 
structure or an extension toward the west from CIS. This 
luminous component is visible in the lower panel of Fig. 6 
but not in the K-map of the top panel. From spectroscopic 
data, K98 argued that this structure may lie at z m 0.3. A 
large void region between the clusters is also apparent in 
the mass map. When decreasing the smoothing scale, the 
core of C1E splits into two maxima that are also visible in 
the higher resolution light map of Fig. 6. 

4.4. Properties of clusters 

The global properties of the three clusters are explored 
using integrated physical quantities enclosed within the 




Fig. 8. Upper panel: Mass reconstruction derived from 
the light distribution emitted by clusters+foreground el- 
liptical galaxies kl ■ This map is almost a reproduction of 
the b) panel of Fig. 6 with one arcminute filtering scale. 
Middle panel: Reconstructed dark matter surface density 
k(6) using the KS93 inversion technique. Lower panel: 
Same reconstruction after 45° rotation of source galaxies. 
No noticeable patterns due to systcmatics are visible. In 
all maps, the shear has been sampled at the observed po- 
sition of background source galaxies, and the maps suffer 
the same (edge+mask)-effects. The Gaussian smoothing 
scale is 1'. Levels are the same as in Fig. 6. The overall 
agreement between panels a) and b) is good. 

radius ro = l/i 70 1 Mpc(= 3'). For each cluster the center 
is set to the X-ray emissivity center. Table 2 summarizes 
the main quantities : the total mass from weak lensing es- 
timates (row 3), the total rest- frame B luminosity emitted 
by the supercluster early- type galaxy sample (row 4), the 
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inferred mass-to-light ratio (row 5), the "spectroscopic" 
velocity dispersion compiled from Dressier & Gunn (1992); 
Fabricant et al. (1994); Carlberg et al. (1996) (row 6), and 
the velocity dispersion derived from a fit of the weak lens- 
ing data to a Singular Isothermal Sphere (SIS) model (row 
7): k{9) = e E /26 , where 9 E = 4tt (f) 2 /? is the Einstein 
radius. Rows (8) and (9) are X-ray ROSAT HRI/IPC and 
ASCA data (Gioia & Luppino 1994; Fabricant et al. 1994; 
Henry 2000) 7 . Row (3) is computed using the densitomet- 
ric statistic : 



C(Mo) = «0' < 9)) - (k(9 < 9' < 9 )) ■ 



( lt {9')) dlnfl'. 



1 - (Oo/OY 



(11) 



M(9) — Y» C ntTr(D d 9) 2 ((9, 9 ) gives a lower bound on the 
mass contained in the cylinder of radius 9. (7*(#)) is the 
average tangential shear. 9o is set to 7 arcmin. In practice 
we used the estimator 



COW 



Var(C) 



(E 4e /^) 2 



(12a) 



(12b) 



with I = {i\9<9i< 9 }. The SIS 9 E value is obtained 
by a x 2 minimisation : 



(13) 



where e t is the tangential ellipticity relative to the cluster 
center. A trivial estimator for 9 E is : 



= 2 



J2i Wi e t,i/9j 

Eiffel 



(14) 



Note however that this estimator is no longer valid when 
k <~ 1. This means that we have to select galaxies far 
enough from centers of clusters. Typically, we set 1' < 
9 < T. 

Finally, the rest-frame B-band luminosity is obtained 
by adding up the luminosities of cluster galaxies with in- 
creasing radius. Systematics due to the selection of super- 
cluster members or to contamination dominate the error 
budget but are small (of order 5% when changing the lim- 
its of Eqs. (1) by 10%). To account for cosmic variance, 
we increased the Poisson noise error by a factor of 1.3, as 
suggested by Longair & Seldner (1979). 

The three clusters differ from one another in terms of 
mass and luminosity. C1N is the most massive and has the 
highest mass-to-light ratio. CIS shows apparent properties 



7 Possible corrections to these values and larger error bars 
may be found in Ellis & Jones (2002); Yee & Ellingson (2003). 
Since the following analysis does not deal with these data, we 
refer to these papers for further information concerning the 
supercluster's X-ray properties. 



of a well relaxed cluster. It is highly concentrated with 
strong lensing features between the two brightest cluster 
galaxies (Mathez et al. 1992) and a rather high X-ray lu- 
minosity. C1E seems more complex: it is the most lumi- 
nous in the i?-band although it is the least massive and 
the least X-ray luminous. Table 2 shows that its kinemat- 
ical velocity dispersion is much higher than what we infer 
from weak lensing. The latter estimate is more typical 
of a cluster mass than the value derived from kinematic 
data. Hence, C1E is likely not relaxed. We attempted to 
describe its bimodal structure (see Fig. 6) by fitting two 
individual isothermal spheres at the location of the lu- 
minosity peaks C1E1 (03 h 06 m 16.5 s , +17° 21' 18") and 
C1E2 (03 h 06 m 19.9 s , +17° 18' 21"). We found cr clE1 = 
312^°° kms" 1 and ct C i E2 = 473l^ kms~ 1 . The fit qual- 
ity is slightly improved, though the quadratic sum of these 
individual velocity dispersions is comparable to the sin- 
gle isothermal sphere fit in table 2. Note that C1E is at 
z = 0.418 which is a rather high radial distance to the 
other clusters. The previous studies of Fabricant et al. 
(1994) and K98 demonstrated that C1E might not be grav- 
itationally bound to the supercluster system. 
The X-ray luminosity presents a better correlation with 
mass than with B-band luminosity. The mass-to-light ra- 
tios are rather different but the mean value within 1 mega- 
parsec is M/L = 249+32 h 70 (M/L) Q . Within 500 kpc 
we found M/L = 23ll™ h 70 (M/L) Q showing that no sig- 
nificant variation with radius is observed. It is worth notic- 
ing that values of M/L for individual clusters have a larger 
scatter. 

We found larger errors than K98 for M, but our estimates 
are not based on smoothed mass maps from randomly 
shuffled catalogs. We directly used galaxy cllipticities in 
Eq. (12). Hence, our error estimates are more conservative 
and do not suffer edge + smoothing effects (+ uncontrolled 
residual correlations) . 



5. Correlation Analysis 

5.1. Linear biasing hypothesis 

The high signal to noise ratios and the good resolutions of 
the light and mass maps are sufficient to explore how light 
and mass correlate and how these quantities evolve as a 
function of angular scale. The statistical properties of the 
relation between dark and luminous matter components 
can then be analyzed from the cross-correlation of the n 
mass map with the K-from-light map shown in panels b) 
and a) of Fig. 8. 

Let us first assume a simple linear relation between the 
luminosity from early-type galaxies (cluster+foreground) 
and the dark matter component. The construction of the 
k e map for the luminosity of early- type galaxies is de- 
tailed in Sect. 3.2. We compute "light" maps again by 
adopting the same scaling relation as in Eq. (3) with a 
starting mass-to-light ratio M/L = 300 h 70 (M/L) . The 
linear biasing hypothesis between the dark matter conver- 
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Table 2. Summary of cluster properties. M is a lower bound on the cluster mass. L E is the rest-frame blue band 
luminosity. M/L is the mass-to-light ratio. These three quantities are calculated inside r = 1 h?Q Mpc. a ve i is the 
kinematic velocity dispersion compiled from spectroscopic data (Dressier & Gunn 1992; Fabricant et al. 1994; Carlberg 
et al. 1996) whereas agis is the velocity dispersion deduced from weak lensing when fitting an isothermal profile for 
the cluster dark matter halo. Note the two distinct values of cr ve i for CIS. The lower value from Carlberg et al. (1996) 
is based on a larger galaxy sample. It is also in better agreement with our estimate, -Lx.boi is the bolometric X rays 
luminosity and Tx the gas temperature (Gioia & Luppino 1994; Fabricant et al. 1994; Henry 2000). 
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gence fields km and ke simply reads : 

k m = Ak e . (15) 

Hence, 300A is the mean mass-to-light ratio. If we assume 
that it is constant with scale and redshift, A is easily con- 
strained by the cross-correlation analysis. 

We compute the two-dimensional and azimuthally av- 
eraged cross-correlations: 

C AB {6) = (k a (#)k b (# + 9)) = (k a k b ). (16) 

We have to subtract the noise contributions to the corre- 
lation functions. Since noise properties of mass and light 
are not correlated, we only have to calculate the noise au- 
tocorrelations : 

(k A K A ) > (kaKa) - (KAKA)nois C - (17) 

Noise autocorrelation as well as error bars are calculated 
by a bootstrap technique. We performed 32 randomiza- 
tions of background galaxy catalogues that mimic the 
noise properties in km as predicted by Eq. (10). We also 
randomly shuffled the shear catalogue calculated with Eq. 
(2) before smoothing and performing the 7-to-K inversion. 
Note also that we discarded the pixels of the convergence 
maps that lie inside masked areas (see Fig. 4). In these 
regions, the lack of background galaxies severely increases 
the noise level. Field boundaries are masked in the same 
way. 

In the following, {k 2 m }, (k e ), and (km^e) refer to the 
mass-mass, light-light and mass-light correlation func- 
tions respectively, (kmke) shows a maximum at zero lag, 
which is significant at the 10-ct confidence level. The cross- 
correlation peak is fairly isotropic and well centered on the 
origin. At zero lag, the normalization parameter of Eq. 
(15) yields M/L = 277± 27h 7 o (M/L)©. 
We thus increased the number of constraints by consider- 
ing the whole correlation function profile over the 7 inner 
arcminutes. The A value is derived by performing a global 



X 2 minimization over the correlation functions, using suf- 
ficiently sparse sampling points to reduce the correlations 
between bins 8 . A satisfies the system: 

(«m)-a 2 (4) 

{k m ke) =\{k e ). 

We found M/L = 2861^ h 70 (M/L) with x7dof ^ 0.88. 
The left panel of Fig. 9 shows the (kmke) and (k 2 e ) cor- 
relation profiles with this mass-to-light ratio normaliza- 
tion. We also observe an excess of light autocorrelation 
at 9 > 15' which is the characteristic distance between 
clusters. Note that this bump is enhanced if we only con- 
sider supercluster early-types and discard the less clus- 
tered foreground contribution. 

So far, we find an excellent matching between the 
(k m ), (km^e) and {k 2 e } correlation functions profiles up 
to ~ 10 arcmin. The linear relation (15) turns out to be a 
good model. As already pointed out by K98, the main 
conclusion is that early-types galaxies trace the mass. 
Oscillating patterns around the light autocorrelation ap- 
pear for r > 8'. G02 as well as Wilson et al. (2001b) found 
similar patterns. They are likely noise artifacts. 

As compared to the results of Sect. 4.4, the correlation 
analysis gives a value for the mass-to-light ratio ~ 15% 
higher than that deduced from integrated quantities in- 
side one megaparsec around clusters. The M/L deduced 
from k maps is insensitive to a constant mass sheet (the so- 
called mass-sheet degeneracy). Therefore, it is necessary 
to subtract the mean luminosity contribution in the cir- 
cular aperture of individual clusters analysis and to only 
consider the excess of luminosity. We find that within 
3 arcmin from the center M/L — 273 ± 47 h-?o (M/L)©. 
Therefore, the agreement with the overall correlation anal- 
ysis is excellent. 

8 1 arcmin is the characteristic length of our spatial smooth- 
ing. We checked that the crossed terms in the covariance matrix 
drop significantly beyond this scale. 
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Fig. 9. Left panel a) : ( k m)) ( k mke) and (n%) correlation functions for a mass-to-light ratio M/L = 286 ± 
36 hjo (M/L)© that fits the correlation functions at scales 9 < 8'. Right panel b) : Same plot with a TIS halo 
model (with truncation radius s* = 150 h?Q kpc). With the same assumption M oc L and a slightly lower % 2 , this 
model also confirms the general conclusion "light traces mass", provided the truncation radius s* < 200 h^Q kpc. Note 
that the bins are correlated. For clarity, in both cases one fifth of the error bars is displayed for (n 2 M ). This coarser 
sampling roughly shows the required spacing for independent bins. 



The agreement with the K98 results, after rescaling to 
a flat A cosmology, is also excellent. Their conclusion that 
early-type galaxies trace the mass faithfully is therefore 
confirmed by our analysis. Nevertheless, the authors ar- 
gued that they saw little evidence for any variation of M/L 
or 'bias' with scale. K98 addressed this issue by performing 
the correlation analysis in the Fourier space by splitting 
the data into a low and a high frequency bin. They found 
an increase of M/ L ratio with increasing wavelength, rang- 
ing from ~ 180 at scales < 2.5 h 7Q Mpc to ~ 280 beyond. 
The physical meaning of this trend is not clear. Variations 
of M/L ratio with scale likely indicate underlying physi- 
cal changes in the relations between mass and light that 
cannot be interpreted from our simple linear scale-free bi- 
asing parameter A. In the following, we investigate some 
models that may explain the M/L variations observed by 
K98. 



5.2. Changing the dark matter halo profile 

In Eq. (3) , we assumed that dark matter halos of individ- 
ual galaxies have a little extension compared to the weak 
lensing filtering scale, so that they can be modeled as point 
masses with mass proportional to the galaxy luminosity. 
In this section we study how a more complex dark matter 
halo profile may change the conclusions of the previous 
section. 

Let us consider a truncated isothermal sphere (TIS) 
(Brainerd, Blandford, & Smail 1996; Schneider & Rix 



1997). The convergence reads 



2r 



1 



Vr 



(19) 



where s is the truncation radius. When s — > oo, b re- 
duces to the Einstein radius of the singular isother- 
mal sphere (SIS). Assuming a L oc a scaling relation 
(Faber & Jackson 1976; Fukugita & Turner 1991) and 
Mtot — 7rS cr i t 6s oc L, we set 



1/2 



L 



1/2 



(20) 



This empirical parameterization is consistent with Wilson 
et al. (2001a) who assumed b oc L 1 / 2 , as well as 
with Hoekstra, Yee, & Gladders (2003) who found b cx 

L 0.60±0.U and s K £0.24+^ leading tQ M K £ 0.84±S;» i 

-, we have to constrain the 



Given that b* = (M/L)^h. 
pair (M/L, s»), or equivalently (A = (M/L)/300, s*). A no 
longer contributes linearly to the he expression because 
of the dependence of s on L. 

The correlation functions are calculated in the same 
way as in Sect. 5.1. However, since s is different from 
one galaxy to another, the resulting correlation function 
(kmKe)tis (resp. (k e )tis) is no longer the convolution of 
(kmKe) (resp. (k\)) by the normalized halo profile (resp. 
normalized halo profile autocorrelation) , making the CPU 
cost much more important. 

Contour plots for (M/L, s») are displayed in the left 
panel of Fig. 10 yielding (M/L) T is = 305±jjj? h 70 (M/L)© 



and s» 



150 



+90 l-1 
-150 "-70 



/i 70 kpc. This value is smaller but still 
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Fig. 10. TIS halo modeling. Left : Contour plot showing constraints of mass-to-light ratio and truncation radius of a 
early-type galaxy. The data are consistent with dark matter halos with truncation radius s* < 200 hjQ kpc. Middle 
: same plot when considering the periphery of clusters only, the tendency is reversed and large halos (s* > 400 h^ Q kpc) 
are favored. Right : same constraints interpreted in the a*, M/L plane. The two cases with/without cluster masking 
are overlaied (solid blue / dashed red contours respectively). The horizontal straight line c* = 220kms _1 is a fiducial 
value for local elliptical galaxies. 



statistically consistent with s* = 264 ± 42 h^Q kpc derived 
by Hoekstra, Yee, & Gladders (2003). However, Hoekstra, 
Yee, & Gladders (2003) used both early and late type field 
galaxies and also relaxed the constraint M cx L, making 
a comparison with our sample difficult. However, because 
we are using a smoothing scale 9 S = 40" = 220 K^q kpc, 
it is only possible to put an upper limit s* < 200 h^Q kpc. 
We therefore cannot rule out that tidal stripping effects 
in dense environments may decrease the galaxy cut-off ra- 
dius, as reported by Natarajan, Kneib, & Smail (2002). 
This point will be discussed in more detail in the next 
sub-section. 

It is also interesting to interpret our results in terms 
of halo velocity dispersion er* = C\J J^* Dd as shown in the 
right panel of Fig. 10. The Results are consistent with 
general values for er» (see e.g. Seljak 2002, and references 
therein). 

5.3. Large scales / Periphery of clusters 

The results derived in Sect. 5.1 and Sect. 5.2 are in 
good agreement with those of Sect. 4.4. They confirm 
that the average mass-to-light ratio of halos is M/L w 
300^-70 (M/L) and that early type galaxies are the pri- 
mary tracers of dark matter on supercluster scales. Their 
contribution may however depend also on the local den- 
sity, and the average value we derived could only reflect a 
biased signature of the mass-to-light ratio dominated by 
the three clusters. One could conclude equally well either 
that early-types trace the mass at all scales with a con- 
stant M/L = 300 hjQ (M/L) or that the signal coming 
from clusters is too strong, and hides more subtle details. 
This would explain why K98 reported an increasing vari- 
ation of M/L ratio with increasing scale using two bins of 
low and high spatial frequencies. 

To clarify this, we calculate the correlation functions 
as above, but we discard the central regions of clusters. 
More precisely, we set to zero the inner 3 arcmin around 



each cluster (circles of Fig. 4) to compute the residual 
correlation produced by the larger scale structures, like 
filaments and voids. When considering the periphery of 
clusters only, the amplitude of correlation functions drops 
by a factor ~ 3 showing that most of the signal comes 
from the clusters. 

The constant M/L ratio model with point-mass-like 
dark matter halo provides a rather bad fit: x 2 /dof ~ 2.1. 
This is significantly worse than for the whole field anal- 
ysis. The best fit yields M/L = 276 ± 3O/i 7o (M/L) 
and is plotted in the left panel of Fig. 11. For a TIS 
halo model the goodness-of-fit is significantly improved 
(x 2 /dof ~ 0.65) when constraining (M/L, s*). However, it 
requires s* > 300 kpc and M /L = 280 ± 40. Note that 
large values of ~ 1.5 hjQ Mpc with larger M/L ~ 400 
are also consistent with the data (sec middle panel of 
Fig. 10). The right panel of Fig. 11 shows such an ex- 
tended halo profile. It is worth noting that these solutions 
appear unphysical and may rather indicate that the input 
model is not well suited. 

The fact that halos are more extended outside the cores 
of clusters is also consistent with the tidal stripping hy- 
pothesis. However, as discussed in the following section, 
this conclusion depends on the input model M cx L and 
on the fact that we assumed that all the mass is associated 
with early-type galaxies. In particular, the contribution 
of late-type galaxies has been neglected again. The small 
amount of mass located in low frequency modes in the the 
outer parts of the supercluster may give an indication that 
these modes are not well traced by early-type galaxies. 

6. Discussion 

The MS0302+17 supercluster mass distribution, de- 
rived from weak lensing analysis of background galax- 
ies, matches the supercluster light distribution of its 
early type galaxies. The correlation between them is very 
strong. More precisely the shape of the light-mass cross- 
correlation profile is in excellent agreement with a simple 
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Fig. 11. Same as Fig. 9 but with the clusters centers. A large truncation radius TIS model (right panel) provides a 
better fit than the point-mass halo model (left panel) which does not fit the data well (% 2 /dof ~ 2.1). Notice how 
compared to Fig. 9, the correlation amplitudes drop when the clusters centers are discarded. 



model where dark matter is directly related to light, as- 
suming a constant mass-to-light ratio. 
We confirm the results of K98, with different data sets and 
a larger field, with different hypotheses to derive lenses 
samples and as well as background sources catalogues, 
and by using an independent PSF correction method. 
Therefore, the strong correlation they found is confirmed 
and strengthened by this work. In particular we confirm 
that M/L ~ 300 ± 40h 70 (M/L) with most matter at- 
tached to the early-type galaxies. A generalization of our 
findings to all supercluster systems is premature but it 
is worth mentioning that Wilson et al. (2001b) and G02 
found similar trends in blank fields and the A901/A902 
supercluster, respectively. 

When we introduce dark matter halos in the form of 
truncated isothermal spheres (TIS) we show that the lin- 
ear relation M oc L is still verified and that dark mat- 
ter halos of early-type galaxies must be rather compact 
(s* < 200 kpc) near the cluster centers that dominate 
the signal. We attempted to mask the clusters centers to 
analyse the remaining signal. Removing these regions be- 
fore doing the correlation analysis gives indications that 
galaxy halos are more extended at the periphery of clus- 
ters than in the inner regions. Such a behavior is also 
consistent with the previous result of K98 who measured a 
different M/L ratio when considering low and high spatial 
frequency modes of K and K-from-light maps. This result, 
which does not have a straightforward physical explana- 
tion, together with our halo analysis, can be interpreted 
in two different ways provided M/L is constant and dark 
matter halos follow a nearly TIS density profile: 

1. either most of the mass is attached to early- type galax- 
ies with M/L = 300 hjo (M/L)© and is distributed into 



halos that are more compact when located closer to 
clusters cores (consistent with the tidal stripping hy- 
pothesis); 

2. or, at the periphery of clusters M/L = constant is 
not completely verified. Late type galaxies (which are 
more abundant at the periphery of clusters) or a more 
diffuse dark matter component that does not follow the 
light from early-types in a simple manner may likely 
be an increasingly important mass component beyond 
the cluster scale. 



Do late type galaxies contribute to the supercluster mass? 

The relation between dark matter and light distribution of 
late type galaxies is difficult to derive from our data only. 
Late type supercluster galaxies cannot be easily extracted 
from our galaxy color-color diagrams, because the B and 
V data are not as deep as the R image. Furthermore, the 
color-color tracks of late type galaxies are broader than 
those of early type galaxies and are therefore much more 
difficult to separate from field galaxies. Nevertheless, we 
find that the relation between dark matter and late type 
galaxies is weaker than and possibly different from that 
of the sample of early-type galaxies. The cross-correlation 
profile can be interpreted as if only a small amount of mat- 
ter is associated to these galaxies. A contribution of low 
frequency modes to the correlation functions is not sur- 
prising since there is compelling evidence that late-types 
are much less clustered and less massive than early-types 
(Budavari et al. 2003). 

From a lcnsing analysis point of view, it is therefore ex- 
pected that the convergence Kiate-typcs is localized in low 
frequency modes, on scales that could be similar to the 
CFH12K angular size. A weak lensing analysis on a single 
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field may not be relevant for probing the mass-light cross 
correlation of late type galaxies on the supercluster scale. 
A much larger field of view is likely needed. Gray et al. 
(2003) used COMBO-17 data in the Abell901/902 system 
and found that late-types are basically located in under- 
dense regions, a result compatible with the well known 
segregation effect. 

On cluster scales the three systems C1N, CIS and C1E 
show clear lensing features (arcs or arclets). Their prop- 
erties show a scatter of mass-to-light ratios with various 
morphological aspects (with indications that CIS and C1N 
seem dynamically relaxed and C1E has ongoing mcrging- 
like events). Overall, the supercluster dynamical state may 
indicate that C1E could not be gravitationally bound to 
the whole system. A large number of redshifts in the field 
would be useful to confirm the global dynamical stage of 
this system. 

K98 possibly detected a filament of dark matter connect- 
ing CIS and C1N. We do not confirm this. We just observe 
an elongated structure in the K-maps which is located 
westward of CIS and is likely due to a casual projection 
effect creating a bridge between CIS and a clump prob- 
ably belonging to the supercluster. Another filamentary 
structure extends toward the West of C1N along the field 
boundary. We indeed also observe a visible counterpart 
in the ke maps. Finally, there is no conclusive evidence 
for any detection of filamentary structure in the field of 
MS0302+17. The detection of K98 may be due to resid- 
ual systematics in the PSF anisotropy correction process. 
Note also that G02 observed a filament connecting A901b 
and A901a but it was not confirmed by an optical coun- 
terpart. Indeed, a detection of dark matter dominated fil- 
aments similar to what is seen in numerical simulations 
remains challenging for such lensing studies. 

Finally, we also observe a large under-dense region 
located between the three clusters. Its angular size is 
about ~ 12'. The depression amplitude is AS « -3 x 
10 12 M /arcmin 2 . 

There is an observational issue that needs clarifica- 
tions. In A901/A902, G02 derived M/L ~ 88 h 70 (M/L) Q 
with correlation analysis whereas they found M/L ~ 
140 (M/L)© in apertures around clusters. Wilson et al. 
(2001b) derived a constant M/L = 210 ± 53 h 70 (M/L) 
for their blank fields sample. These values are significantly 
different from K98 and this work. The reason for this dis- 
crepancy is not clear. The large scatter in M/L found by 
G02 from one cluster to another may be intrinsic if one as- 
sumes that each cluster is in a different dynamical stage. 
G02 investigated whether the large scatter could be in- 
terpreted as a natural scatter in the mass/light relation. 
Using the Dekel & Lahav (1999) formalism, they claimed 
they measured a marginal nonzero stochastic component 
in the Abell901/902 system. In the case of MS0302+17, 
we are unable to measure such a positive stochastic term 
in the correlation function profiles. The subcomponents 
of the A90 1/902 supercluster are physically closer than 
those found in MS0302+17. The average projected physi- 
cal separation of the former is of order 2 hj~ Mpc whereas 



the MS0302+17 clusters are separated by ~ 5 hj~Q Mpc 
showing that possible interactions and dynamical stages 
are different from one supercluster to another. 

7. Conclusions 

We have analyzed the weak lensing signal caused by the 
supercluster of galaxies MS0302+17 and connected it to 
its optical properties. Using a BVR photometric dataset 
from CFH12K images, we identified the early type mem- 
bers of the supercluster. The R band image was also used 
to measure the coherent gravitational shear produced by 
massive structures of the supercluster and by foreground 
contaminating field objects. 

When considered individually, each cluster has an av- 
erage rest frame B band mass-to-light ratio M/L = 
249l*2 ^70 (M/L)©- The Eastern cluster does not show a 
well relaxed structure. It can be viewed more likely as a 
two component cluster system with ongoing gravitational 
interaction. This may explain the rather poor agreement 
between lensing and kinematic estimates of the velocity 
dispersion. It also supports the previous conclusions of 
Fabricant et al. (1994) and K98 that C1E may not be grav- 
itationally bound to the system made of the other two C1N 
and CIS clusters. 

The mass (or convergence) map shows an excellent agree- 
ment with that derived from the distribution of early-type 
galaxies. Besides the well detected main clusters, one can 
observe a large underdense region between them. We were 
unable to confirm the existence of a filament joining CIS 
and C1N as claimed by K98. 

We performed a correlation analysis between "light" and 
mass aiming at probing whether the linear relation Mai 
(or more precisely k oc Kf rom light) is consistent with the 
data at hand. The results based on mass-mass, mass-light 
and light-light correlation functions are robust enough to 
make conclusive statements on the average mass-to-light 
ratio. We found that M/L ~ 300 h m (M/L) . In other 
words, all the mass detected from weak lensing analysis is 
faithfully traced by the luminosity distribution of early- 
type galaxies. 

Our conclusions are in excellent agreement with those 
of K98. They only depend slightly on the unknown dis- 
tribution of late-type galaxies, since their contribution is 
found to be small. However, when focusing on early types, 
we were able to put constraints on the density profile of 
galaxy halos. Despite the rather large spatial smoothing, 
we conclude that halo truncation radii s* < 200 hj^ kpc 
for an L* galaxy. We also found evidence for a relax- 
ation of this constraint at the periphery of clusters: s* > 
300 /if kpc. However, this latter result relies on the fact 
that late-type galaxies are neglected. Such an hypothesis 
may not be so evident at large distance from the centers 
of clusters. 

Further investigations of the MS0302+17 supercluster 
of galaxies may require more photometry (in different op- 
tical/NIR bands) in order to identify late- type galaxies in 
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the supercluster and compare their distribution and phys- 
ical properties with the early-type sample. 
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